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This paper presents a method for simultaneous optimal structural and 
control design of large flexible space structures (LFSS) to reduce vibration 
generated by disturbances. Desired natural frequencies and damping ratios 
for the closed-loop system are achieved by using a combination of linear 
quadratic regulator (LQR) synthesis and numerical optimization techniques. 
The state and control weighting matrices (Q and R) are expressed in terms of 
structural parameters such as mass and stiffness. The design parameters are 
selected by numerical optimization so as to minimize the weight of 
the structure and to achieve the desired closed-loop eigenvalues. An 
illustrative example of the design of a two bar truss is presented. 


INTRODUCTION 

Large structural systems in general and large space structures in 
particular present new challenges to the structural dynamicist and the 
control engineer as well. Indeed, such large systems may exhibit well over 
a thousand vibrational modes usually closely spaced and with little, if any, 
damping. Some form of active control is likely to be necessary in order to 
meet exacting stability and pointing requirements. In fact, structural 
requirements (primarily low mass) increase the need for active control. 
Some optimal trade off between structural and control criteria has to be 
achieved. 

Until recently, the design of control systems for large structural 
systems was a two-step procedures: first the structure was designed based on 
structural criteria (primarily total weight); then in a second step a 
control system (satisfying some desired control objectives) was designed for 
the structure obtained in the first step. Inasmuch as a low weight (and 
thus low stiffness) structure will require high control energy, the design 
objectives of the two steps are to some extent contradictory so that an 
optimal control design for an optimally designed structure will not in 
general result in an overall control-structure optimal design. Both designs 
need to be carried out simultaneously. 
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LITERATURE REVIEW 


The optimal structural and control design of large flexible space 
structures was recently investigated by several researchers. Venkayya and 
Tischler [1-2] have suggested that the performance index (PI) in optimal 
control of structural systems be a measure of the system total mechanical 
energy. By appropriately choosing the state and control weighting matrices, 
the PI can be expressed as the (weighted) sum of the kinetic, strain and 
potential (including control) energies. Knot and Venkayya C 3~ ^ 1 tackled the 
structural and control optimization problem by minimizing the weight of the 
structure with constraints on structural frequencies and the minimum 
Frobenious norm of the gain matrix. This process has to be carried out in 
an iterative fashion. 

Becus and Lui [5] have proposed a general method to choose state and 
control weighting matrices in optimal control design so as to satisfy 
desired closed-loop eigenvalues. This was further extended by Becus and 
Sonmez [6] to allow for eigenvector assignment. In this paper we combine 
both ideas in order to obtain a method to carry out simultaneous optimal 
structural and control design. 

Desired dynamic structural requirements (natural frequencies and 
damping ratios for example) can be expressed both in terms of desired 
closed-loop ei genstructur e (eigenvalues and/or eigenvectors) and structural 
parameters (mass and stiffness for example). Using a PI of the form 
suggested in [ 1 J , the elements of the state and control weighting matrices 
(Q and R respectively) are also expressed in terms of structural parameters. 
Thus, when choosing the Q and R matrices (using the method of [5-6]) to 
satisfy a desired closed-loop ei genstructur e (i.e. dynamic structural 
requir ements ) , one in fact chooses new strucutral parameters and therefore 
carries out a simultaneous optimal control structure design. 

In this paper a new design algorithm is developed so that a minimum 
weight structure with desired damping and natural frequency of the closed- 
loop system can be obtained. We compare the results with [33 in the last 
section. 


SIMULTANEOUS STRUCTURAL AND CONTROL OPTIMIZATION 

Consider a controlled structural dynamic system described by the 
discrete (finite element) model 

Mr + Kr = Du O ) 

where r is a vector of n physical displacements and the number of control 
inputs (forces) u is m. M, K and D are the mass, stiffness and applied load 
distribution matrices of appropriate dimensions respectively. Assume that M 
and K are positive definite. 

The state space representation of Eq. (1) can be written as 


x * Ax + Bu 


(2) 


where 


x 



(3) 


226 



A ( p ) , 
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B(p) , 


where p is a vector of structural parameters of dimension 1. 
The optimal steady-state control is a linear state feedback 


(M) 


(5) 


u - -Gx . (6) 

The state feedback gain matrix G is obtained from LQR synthesis and the 
closed-loop system is given by 


x - (A - BG)x 


(7) 


LQR synthesis determines a control u which minimizes the quadratic 
performance index [1] 

PI * fo [o m r T Mr + Q k r T Kr + 0 r u T D T K _1 Du] dt (8) 


or in the state space coordinates 



PI * 

/„ [x T Qx + u T Ru] dt 

(9) 

where 

Q - 

0 M 1 0 

m ' 

* Q(p) 

(10) 

0 I V 

and 

R - [ 

i dV'd] ■ 

r J 

R(p) 

(11) 


for positive scaling parameters © m , o k and 0^ . In Eq. (8), PI is the 

absolute weighted sum of the kinetic, strain and potential energies. 

The relationship between characteristic Dolynomial of the optimal 
system and weighting matrices is obtained as follows [7] 


or 

where H(s) 


det(sl - Z) 


si - A 

br -1 b t 

Q | 

1 Si + A T 


* c (s)<l>c(- 3 > ■ <l>o ( 3 ) <f> o ( ~ 3 ) de t [ I +R~ 1 H T ( -s ) QH ( s ) ] 
is the open-loop transfer function matrix 


( 12 ) 

( 13 ) 
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H(s) - (si - A ) “ 1 B , (^) 

Z is the canonical system matrix, <}> c (s) anc ^ 4>o(s) are the closed-loop and 

open-loop characteristic polynomials respectively. 

For a given desired closed-loop pole s * which is not an open-loop 

pole, the determinant in the right-hand side of Eq. (13) must equal zero 
when the weighting matrices Q and R take values which yield the desired 
closed-loop eigenvalues. In order to use numerical optimization techniques 
to solve Eq. (13) for Q and R, we, as in Ref. [8], set the objective 
function as 


obj = det[l + R _ 1 H T (-s d )QH(s d )] - 0 . (15) 

The desired characteristic equation corresponding to Eq. (15) is 


* 11 4 

J *1 


(s - s Ms 


S d. } 

J 


(16) 


where s J is the j-th desired closed-loop eigenvalue, 
j 

Q and R are determined by equating coefficients of the terms involving equal 
powers of s in Eqs. (15) and (16). This yields 


f,(p) = 0 


(17) 


f k (p) - 0 

where k is the number of equality constraints which involve equal powers of 
s in Eqs. (15) and (16). 

The objective in structural and control optimization is to make the 
selection of design parameters so that the structure weight is a minimum and 
the specified closed-loop eigenvalues are satisfied. The optimization 
problem can be stated as 

Minimize the weight W * W(p) 

subject to Eq* (17) (18) 

and P g > P g , s - 1 , • • • , 1 , 

where P denote minimum allowable values of the structural design 
s 

parameters . 
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ILLUSTRATIVE EXAMPLE 


In order to illustrate the feasibility of the above algorithm, the 
structural two bar truss model shown in Fig. 1 was considered as a simple 
design example. For the geometry shown, the dynamical equations of motion 
(Eq. (1)) are 


2 

0 




(A| +A 2 ) 

2(A,-A 2 ) 




cos( 6) 





+ k, 
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sin( 0) 


thus M 
and 
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2 0 
0 2 
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(A, + A 2 ) 
2( Aj - A 2 ) 


2 ( A , - A 2 ) 
4(A, + A 2 ) 


( 20 ) 


( 21 ) 


are the optimal mass and stiffness of the structure respectively. In Eq. 
(19), A, and A 2 are the cros3-3ectional areas of the bars and k, - E/(5L) is 

a stiffness coefficient, E representing the elastic modulus of the bars and 
L the length of the members. A control force u is located at the vertex 
with 9 being the angle between its line of action and the horizontal. r, 

and r 2 are the horizontal and vertical displacements of the vertex 


respectively. 

The dimensions of the structure were given in unspecified consistent 
units. The elastic modulus of the members was assumed to be 1 and the 
density p of the structural material was assumed to be 0.001. A 
nonstructural mass of 2 units was attached at node 2 and the structural mass 
of tne members was ignored for simplicity (thus the mass matrix of Eq. 
(20)). The actuator and sensor were located in element 1 connecting node 1 
and 2. The minimum cross-sectional area was set equal to 10 units for both 
members. 

Once the choice of the material is fixed, the design variables are the 
cross-sectional areas of the members Ai and A 2 , the scaling parameters 0 m , 

and the angle 0 of the applied load with respect to the 


0, and 0 , 
k r ’ 


horizontal. The optimal closed-loop eigenvalues are specified as s^ - - 

d j 

0.0228 + 1 . 1 7j and s^ * -0.361 ± 4.8lj. Arbitrary lower and upper values 


of 0 were 3et at 30° and 60° respectively. 

TM 

Analytical and numerical computations were carried out using MACSYMA ‘ 
[93 for symbolic algebraic manipulations, MATLAB [10] for matrix 
computations and LQR synthesis, and GRG2 [111 for numerical optimization. 
The numerical results for several representative optimal designs are listed 
in Table 1. A discussion of these results appears in the next section. 
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Table 1 . Optimal Two Bar Truss Designs 





































DISCUSSION OF RESULTS 


For this simple example tiie six design variables were not independent. 

The scaling parameters 0 and 0, appeared in the constraint equations only 

m K 

as the combination k,0 /0, . This combination was then used as one of five 

m k 

independent design variables. To obtain the values of Table 1, 0 m was 

arbitrarily set equal to 1 then 0^ was evaluated by multiplying the value 

obtained by numerical optimization by k, » 0.0089. 

Since there are five independent design variables and only four 
equality constraints, there are many solutions to the optimization problem. 
In order to obtain a unique solution one could arbitrarily fix the value of 
one of the five independent design variables or equivalently introduce an 
additional constraint. 

Of all designs presented in Table 1, Design 6 is the best since it 
leads to the lowest value for the weight. This "optimal" design leads to a 
weight of 2.84 which is less than half of the best design of Ref. [3] 
(6.417). 

A closer examination of Table 1 leads to some interesting observations. 
Designs 3 through 6 have weights which are inversely proportional to © k » In 

fact the product O^xWeignt is nearly constant for these four designs and 

equal to 24.45. In addition it can be seen that for these four designs the 
ratio Ai/A, is nearly constant and equal to 10. It is conjectured that many 
other designs could be obtained by choosing areas satisfying this 
relationship and calculating the corresponding 0^ while keeping the other 

design variables constant. 

Design 1 is representative of several designs for which the ratio A,/A 7 
is nearly constant and equal to 0.1 while Design 2 leads to an angle less 
than the upper bound value of 60° . For all designs obtained the product 
0 k *Weight was nearly constant and equal to 24.45. 

Finally it must be noted that as more weight is given to the control 
effort the achieved closed-loop eigenvalues are closer to the desired 
eigenvalues. As more weight is being given to the strain energy cost the 
total weight decreases. 


CONCLUSION 

An algorithm for simultaneous structural and control optimization 
design of a minimum weight structure with desired closed-loop eigenvalues 
was proposed. It has been shown that structural and control designs can be 
obtained by LQR assignment. The design parameters were appropriately 
selected by numerical optimization so as to minimize the weight of the 
structure and to achieve desired natural frequencies and damping ratios. 
The feasibility of the algorithm was demonstrated by applying it to a simple 
example. Further work is needed to investigate the application of the 
algorithm to large-order systems. 
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